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Many studies have confirmed that cracks in brittle materials branch when the crack speed exceeds 
a certain threshold, Vc, but the value of that threshold is not understood. Almost all theoretical 
calculations overestimate Vc by factors of two or more. We show that thermal noise can reduce the 
threshold by a substantial amount, and we propose that this effect can account for the discrepancy. 



The formation of cracks in the fracture of solids leads to 
extremely intricate pattern formation. One of the most 
interesting phenomena in this field is the mirror-mist- 
hackle transition in brittle fracture, for example, that of 
glass . This occurs when a crack accelerates and passes 
through a critical velocity, Vc above which it becomes 
unstable against branching Below Vc the surface 

created by the crack is smooth (mirror-like) while above 
it is rough. However, the value of Vc is not theoretically 
understood: typically, theoretical estimates are too large 
[^,0. In this paper we will show that including thermal 
noise in the theory can lower the threshold substantially. 

For many conditions, once started, cracks accelerate 
until their velocity is of the order of the free surface wave 
(Rayleigh) velocity cr. In the case of a Mode I crack, the 
crack cannot move faster than the sides can move apart: 
that is, faster than cr. Thus cr is an upper bound on the 
crack speed However, in practice, cracks travel slower 
than cr. This has been associated with an instability 
against branching since the seminal work of Yoffe |^ . She 
showed analytically that a straight crack is intrinsically 
unstable at speeds above w 0.66cfl;. The idea is roughly 
as follows: at low speeds there is stress concentration in 
front of a crack leading to further motion. However, at 
high speeds the stress pattern becomes more isotropic, 
and above a velocity threshold the maximum stresses are 
at a finite angle to the direction of the tip; this may be 
assumed to give rise to branching. The effect is due to 
the fact that the stresses are transmitted by sound waves: 
if the speed of the crack is too large the crack 'catches 
up' with the stresses. Other views of the branching in- 
stability have been formulated which differ in detail, 
but not in principle, from Yoffe's work. Simulations of 
2-dimensional spring models also give a branching 
instability around 0.66c/{, and the model we use here 
has the same property. However, the experiments always 
show branching at still lower speeds Vc ~ 0.36cr in the 
case of PMMA, a nearly ideal brittle material; similar 
observations have been made in ordinary glass It 
appears that something is causing branching before the 
Yoffe mechanism takes hold. Here we introduce a phys- 
ical effect that has been left out of previous treatments, 
thermal noise, and show that it could resolve the discrep- 



ancy. 

It is obvious that large enough thermal noise could 
cause a crack to wander. The essential question is 
whether reasonable values for the parameters give rise to 
a substantial effect. A priori this seems unlikely: usually, 
bonds break due to thermal effects only near the melting 
point, which is not the situation at hand. However, this 
is not correct because the nature of the branching is that 
that even below the threshold the stress pattern near the 
crack tip is approaching instability, so that a relatively 
small thermal fluctuation could break the 'wrong' bond, 
and give rise to branching. We will demonstrate this ef- 
fect in a simulation where we attempt to fit all of the 
parameters to experimental values. 

To account for thermal noise, we need to have a cou- 
pling between the long-wavelength elastic modes, which 
we describe by a spring model, and a 'heat bath', i.e. the 
rest of the degrees of freedom of the solid. This arises 
from a damping or viscosity term. In the macroscopic 
theory of elasticity damping can be included in the elas- 
tic equations of motion by introducing Kelvin viscosity 
terms. With these terms the equations of motion for the 
spring model we use are: 

mdttUr = ^ Kh [h ■ (Ur - Ur')] 
r' 

-\-^'nh[h- {dtUr - dtUr' )] (1) 
r' 

where Up is the displacement of the node at r , K and 
Tj are the spring constants and viscosity respectively, n 
is the unit vector from node r to r', and the r"s are the 
nearest neighbors of r. The viscous terms are the simplest 
form consistent with the symmetry of the situation [p^ . 
We have previously shown that if j] is of order unity this 
term is sufficient to allow a crack to propagate without 
branching at a terminal velocity less than Vc. Without 
viscosity any crack, once initiated in in fixed grip loading 
conditions (which is the case we consider) , accelerates to 
the branching threshold |jll|,|l2| . 

The value of t] is central to our theory. We get it by 
fitting to the experimental values for the attenuation of 
longitudinal sound jl^]. It is easy to show from Eq. (|l|) 
that a{uj), the inverse of the attentuation length for lon- 
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gitudinal sound of frequency ui, is given by {ri/2KcL)uj'^ 
in our spring model. We define a characteristic frequency, 
— cl/cl for any elastic system, where cl the speed of 
longitudinal sound, and a is the length unit, which we 
take to be the radius of the crack tip. (It is well known 
that any fracture model must introduce some microscopic 
length unit.) Then we put aa = Q{uj/lJo)'^ ■ The dimen- 
sionless number Q is a measure of the size of the damping. 
For the spring model, in units where K = m = a = 1 
we find that Q = O.5377. For Si02 the sound attenua- 
tion has also goes like over a large frequency range 
iQ. Using the measured data in Ref. for a and cl 
and using a « 0.5 nm, a molecular dimension, and we 
find Q 0.4. For the case of PMMA the attenuation 
does not have the same frequency dependence and 
we cannot unambiguously determine 77 . In the re- 

mainder of this paper we will mostly discuss Si02. For 
this case we are justified taking 77 to be of order unity. In 
all of our simulations we take rj — 2. 

In order to simulate fracture we take Eq. (|l|) as the 
dynamics of a snapping spring model: If any bond is 
stretched more than uth = 0.1a we break it. Strips of 
lattice along the x direction are stressed by holding two 
edges parallel to the y-axis at fixed displacement. A, and 
a crack is started parallel to the x-axis. We artificially cut 
the springs at an effective crack velocity of 0.05 [ pl]JT2| 
at the center of the width along the length. We solve the 
equations of motion by using a Verlet algorithm. As soon 
as velocity exceeds 0.1 the artificial scissoring is stopped. 
The results are unaffected by this procedure. The width 
(y-direction) along which the strain is applied is up to 
100a and the length (x-direction) , along which the crack 
propagates is up to 400a. For more details of our proce- 
dure see [|ll[l2|. 

We first show that the Yoffe instability, or something 
very much like it, is in fact present in our model. We 
create cracks at various terminal velocities by varying A, 
and measure the ratio of the strain in the 'right' bond 
(corresponding to a straight crack) and the 'wrong' one. 
This is shown in Fig. (|l|). The ratio goes through unity 
and the crack branches at « 0.7cr. 

Now we add thermal noise by putting a random force 
term, 9{r,t), on the right hand side of Eq. (|l|). The 
components of the force are independent gaussian noises 
at each site whose amplitude is chosen to be consistent 
with the fluctuation-dissipation theorem: 



2kBTr] 
dt 



St.t'Sr.r'Si. 



(2) 



where i, j label the components. The dependence on 1 /dt 
is the standard way to normalize a discrete (5-function. 

Strictly speaking the noise should have nearest neigh- 
bor correlations in order to be consistent with our form of 
77. Put another way, short wavelength modes have more 
damping than long wavelengths as mentioned above, and 
we should drive them more. Since bonds are broken 



by short wavelength fluctuations, we are underestimat- 
ing the effect of noise by neglecting the correlations. We 
have done a few runs with the completely consistent noise 
to verify this. The difference is not very large. 

We will need to know the energy unit for our model, 
Ka^, in order to evaluate the noise strength. This we 
get by fitting Ku'^yj2a^ to the fracture energy, 8J/m^ 
for Si02 |jl|. Taken literally, this would give each bond 
about 20 cV, which is too large by a factor of 10 or so. 
In our units room temperature comes out to be of order 
10-5 - 10-6. We put ksT = 3 ■ IQ-^. By taking a large 
energy unit we are once more underestimating the noise 
strength. 

In the vicinity of the tip there are large local heat- 
ing effects in some materials including PMMA |l^. To 
account for this we start the lattice with an ambient tem- 
perature (room temperature) and allow it to locally heat 
(and diffuse heat) from the work of the viscous forces. In 
that case the temperature becomes a local variable T(r). 
In order to find r(r) we set: 



dtT{Y) 
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Here k is a thermal conductivity, C is the heat capacity 
associated with the local heat bath, and is the viscous 
force on node r (the second term on the right-hand side 
of Eq. (|l|) . If C is very large the lattice remains at the 
ambient temperature. In our simulations we adjust C 
and K so that the vicinity of the tip heats up to about 
three times room temperature. 

The simulation precedure is similar to that in the noise- 
free case outlined above. We use a version of the Ver- 
let algorithm modified to account for the random forces 
6i{r,t). An initial crack of length 50a is created and the 
lattice is allowed to thermalize. After the temperature 
has settled down, the lattice is strained to the desired 
value. 

Our results are as follows. Since most of the detailed 
observations in the literature are for PMMA, we will com- 
pare to the qualitative features seen in that material, even 
though the quantitative fitting that we have done is for 
glass. Fig. (||a) shows the effect of thermal noise on the 
crack velocity. Noise leads to fluctuations in the crack 
velocity, an effect that is experimentally observed [^,||. 
The large change in the velocity in the figure marks the 
change in the row of propagation. As is seen from the 
figure, there is small rise in the velocity before it satu- 
rates. This is also observed when the fracture energy is 
low In Fig. (||b) we show two cracks, one without 
noise and one with, in otherwise identical conditions. 

With noise there is a finite probability of wandering 
for any v. Thus we must define the critical velocity in 
terms of the roughness of the crack. We call v^, the veloc- 
ity above which the roughness is more than 10% of the 
width. Beyond this point the fluctuations grow quickly 
and saturate due to the finite size of the sample. This 
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observation is consistent with the experiments ||]. Fig. 

shows the effect of thermal noise on Vc- Here the rms 
deviations are measured from the change of direction and 
are averaged over 5 runs for each point. The critical ve- 
locity for crack wandering is lower than in the absence of 
noise and is around 0.45cfl;. 

The effect of crack wandering from one plane to an- 
other causes the rough interface. The experiments ||] 
show that at higher velocities the roughness of the 
cracked surface is larger. Similar effects occur for 
simulations of crack propagation in silicon. The nature 
of the wandering instability on our simulations is inter- 
esting. Just above Vc the crack does not change its di- 
rection of propagation. In our simulations this occurs for 
velocities higher than 0.7cfl. At lower v the crack wan- 
ders by jumping rows to planes adjacent to the initial 
one. That is, as the velocity increases in the presence of 
noise, the stress field component ayy fluctuates and be- 
fore a^x becomes comparable to ayy the crack ceases be 
confined to the initial plane. In Fig. we plot relative 
displacements in the y direction for two cases. Fig. (^a) 
shows these displacements just before the crack begins 
to propagate. Fig. (^) is same plot but just before the 
first wandering is encountered. In the second case the 
relative displacement in y -direction is higher in the row 
close to the middle one in which crack is propagating. 
Recent experiments Q clearly show similar observations. 
There are two velocity thresholds, Vci and Vc2 such that 
for velocities Vci < v < Vc2 surface roughens with almost 
negligible branching while above Vc2 it branches. The 
value of Vc2 is around 0.67cr, (the Yoffe velocity) where 
becomes comparable to Gyy. 

In Fig. (j^) we show as a function of the heat capac- 
ity of the heat bath. When C is large the tip of the crack 
remains at the ambient temperature. We see that about 
half of the reduction of Vc is due to tip heating for our 
parameters. In this paper we have used a naive model 
to treat fracture in Si02. In fact, Griffith used exactly 
the same sort of considerations in his pioneering work 
on glass Q . We think that the qualitative results of this 
study will survive a more detailed look at the theory, and 
that the effect of thermal fluctuations on is potentially 
quite large. The real test for our ideas must come from 
experiment. For example we could imagine looking at the 
branching threshold systematically as a function of tem- 
perature, or, in materials like PMMA where the polymer 
chains vary in length, as a function of heat capacity. 
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FIG. 1. The branching instability as a function of v. The 
strain in the row corresponding to forward growth is compared 
to that in the next row. 
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FIG. 2. Effect of noise on crack propagation, (a.) Dashed 
line is the velocity for a 50 x 400 lattice for T=0°K, r?=2.0, 
specific heat C — 100, k=4.0, and stretch per unit width 
A/VF=0.033 where W is the lattice width. The solid line is 
the velocity for the same parameters with T = 3.0 ■ 10~^ (b.) 
Crack pattern with and without noise in the same loading 
conditions 
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FIG. 3. Rms deviations in the crack path as a function of 
velocity for t] = 2.0, C = ICQ, k = 4.0, T = 3.0 • 10~^ A is 
varied to obtain the appropriate velocity. 
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FIG. 4. (a) Relative displacements in the columns (along 
y direction) containing the crack tip and one column ahead 
just before the crack begins to propagate, (b) Relative dis- 
placements in the columns containing the tip just before first 
wandering is encountered. The parameters are the same as in 
the previous figures. 
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FIG. 5. Critical velocity as a function of heat capacity. For 
large heat capacity local heating in the vicinity of the tip is 
suppressed. 
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